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Formation and decay of rod-like and slab-like nuclei, as may be encountered in the outer 
parts of neutron stars, are examined within a zero-temperature stability analysis with respect 
to perturbations inducing fission and proton clustering. The nonspherical nuclei are found 
to be stable against these perturbations, implying the possibility that the size of the stellar 
region containing such nuclei exceeds the equilibrium prediction. 

Nuclear matter exhibits the coexistence of a liquid phase with a gas phase at 
00 ' subnuclear densities, due to the attractive force responsible for nuclear binding. 

This feature has been revealed in the vicinity of the critical point of the liquid- 
gas transition through the detection of multifragmentation occurring in heavy-ion 
collision experiments performed at intermediate energies.© Such a coexistence is 
believed to be present in the outer part of a neutron star.!! Since the gravitational 
coupling is far smaller than the Coulomb coupling, the star's gravitational stability 
requires nuclear matter to be neutralized by electrons. The resulting electrically 
neutral matter is usually referred to as neutron star matter. At sufficiently low 
temperatures, relevant to neutron star interiors, long-range Coulomb interactions 
force the system to separate periodically into liquid and gas segments of macroscopic 
size, adding a crystalline property to the liquid-gas coexistence. This liquid-gas 
mixed phase, having two conserved charges, the baryon number and the electric 
charge, is manifested for a finite range of pressures. The presence of electrons controls 
the proton fraction in the liquid and gas phases: As the electron chemical potential 
increases with increasing pressure, a system in (5 equilibrium has the liquid phase 
neutron-enriched and leaves the gas phase almost free of protons. 

At zero temperature, it is believed that the energetically favourable configura- 
tion of the mixed phase possesses interesting spatial structures; %%> the liquid (gas) 
phase is divided into periodically arranged parts of roughly spherical, rod-like or 
slab-like shape, embedded in the gas (liquid) phase and in a roughly uniform elec- 
tron gas. Hereafter, we refer to these parts composed of the liquid and gas phases 
as 'nuclei' and 'bubbles', respectively. Recent calculations, Derformed within the 
Wigner-Seitz approximation using specific nuclear models, tP'EP'LP indicate that at a 
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density of about 10 g cm -3 , rather small compared with the saturation density of 
symmetric nuclear matter, p s ~ 2.7 x 10 14 g cm -3 , the spatial structure changes from 
a bcc Coulomb lattice of roughly spherical nuclei to a two-dimensional triangular lat- 
tice of rod-like nuclei. With the density increased further, it is transformed into a 
layered structure composed of slab- like nuclei and bubbles. Next, a two-dimensional 
triangular lattice of rod-like bubbles and a bcc Coulomb lattice of roughly spherical 
bubbles appear in turn. Finally, at a density of about p s /2, the system dissolves into 
uniform nuclear matter. These changes of the spatial structure, accompanied by a 
reduction of the total surface area, are governed by the competition between the 
electrostatic and surface energies. Since slabs and rods look something like lasagna 
and spaghetti, the phases with positional order of one and two dimensions are often 
referred to as nuclear "pasta" . 

In a separate paper,™ we more extensively calculated the corresponding equi- 
librium phase diagrams for zero-temperature neutron star matter at subnuclear den- 
sities. In this calculation, we used a compressible liquid-drop model developed by 
Baym, Bethe and PethicklH' (hereafter denoted BBP), into which we incorporated 
uncertainties in the surface tension E SUI f and in the proton chemical potential p p 
in pure neutron matter. With an increase of E suv {, accompanied by an increase in 
the sum of the electrostatic and surface energies, the density p m at which the system 
becomes uniform decreases. There is also a tendency for p m to decrease with decreas- 
ing pp \ This is because —pp°^ represents the degree to which the gas phase favours 
the presence of protons in itself. It was shown that while the phases with rod-like 
bubbles and with spherical bubbles can occur only for unrealistically small E snr f , the 
phases with rod-like nuclei and with slab-like nuclei survive almost independently of 

For these two nuclear phases, there is a direction in which the system is trans- 
lationally invariant. As noted by Pethick and Potekhin,0) this situation is geomet- 
rically similar to a liquid crystal rather than to a rigid solid. The elastic properties 
of the nuclear rods and slabs can thus be described by elastic constants used for 
the corresponding liquid-crystal phases, i.e., columnar phases and smectic A phases, 
respectively. Pethick and Potekhin expressed these constants in terms of the elec- 
trostatic and surface energies. 

We predicted in Ref. ||) that the phases with rod-like nuclei and with slab- 
like nuclei are energetically favoured in the density regime just below p s and at 
zero temperature, and it is thus important to consider how these nuclei are formed. 
Such formation requires simultaneous migration of an infinite number of nucleons, 
in contrast to the case of ordinary chemical reactions. This prevents the nuclear 
system from crossing the energy barrier formed between the initial and final states 
in the configuration space via quantum tunnelling. Instead, "pasta" formation can 
be driven by instabilities with respect to fluctuations around the initial state. 

In this paper, we examine the kinds of instabilities that are involved in the 
formation and decay of rod-like and slab-like nuclei at zero temperature. As such, we 
first note an instability with respect to quadrupolar deformation of spherical nuc 
as originally investigated in the context of nuclear fission by Bohr and Wheeler. 
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Pethick and Ravenhall& predicted that this instability occurs in the bcc Coulomb 
lattice and creates elongated nuclear rods when the volume fraction of the liquid 
phase reaches about 1/8 in the course of compression. We also consider an instability 
with respect to proton clustering in uniform matter, as considered by BBP.I) This 
clustering is induced by the isospin symmetry energy. At an instability point that the 
system reaches during decompression, the gain due to this energy compensates for 
the gradient and Coulomb energies produced by the resulting inhomogeneities, and 
hence a phase with nuclei of some form appears. Pethick and RavenhallH^ showed 
that the critical densities for instabilities with respect to quadrupolar deformation 
and proton clustering are close to the corresponding equilibrium transition points. 

We extend these considerations to other changes in nuclear shapes. The possible 
instability with respect to proton clustering in planar and cylindrical nuclei tends 
to divide each nucleus into nuclei of lower dimension, while the possible fission- 
like instability of slab-like and rod-like nuclei tends to lead to the formation of 
uniform matter and slab-like nuclei, respectively. Using a typical nuclear model, 
we find that planar and cylindrical nuclei are stable with respect to deformation- 
induced fission and proton clustering. This, together with the finding that these 
nuclei are thermodynamically stable and do not exhibit proton drip, suggests the 
possibility that they persist beyond the equilibrium transition points, e.g., up (down) 
to a critical point at which the stable uniform (bcc) phase nucleates via quantum 
tunnelling of the energy barrier in the configuration space. (Note that a critical 
droplet of the stable phase considered here, which, after forming, develops into bulk 
material, has a finite size, in contrast to the critical droplets of the "lasagna" and 
"spaghetti" phases.) Finally, the implications of such persistency for neutron star 
structure are discussed. 

We begin_by recalling the mechanism of nuclear fission investigated by Bohr 
and Wheeler. t3 They regarded a nucleus as a spherical liquid drop of radius R and 
total charge Ze in which neutrons and protons are distributed uniformly, and then 
calculated the change in the sum of the surface and electrostatic energies induced 
by various kinds of deformations. Since we are interested in the limiting case in 
which the fission barrier vanishes, it is sufficient to consider a slight quadrupolar 
deformation of a spherical drop, characterized by the distance from the center of the 
drop to an arbitrary point on the surface with polar angle 8: 

X{6) =#[1+00 + 0^2 (cos #) + •••] . (1) 

Here, ao is the fractional change in the mean radius of the surface, and a2-P2(cos#) 
represents the degree of the quadrupolar deformation. Nuclear saturation leads us 
to assume the volume of the drop to be invariant under the deformation (Q); we 
may thus write ao = — (l/5)a|. From the condition that the resulting change in the 
sum of the surface and electrostatic energies is zero, we can derive the well-known 
fission-instability relation 

4 0) = 2E f ] ■ (2) 

Here E$ = (3/5)Z 2 e 2 /R is the Coulomb self energy, and £ s (0) = 4irE surf R 2 is the 
surface energy with the surface tension E SUI f. 
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In neutron star matter at zero temperature, these drops form a bcc lattice em- 
bedded in a roughly uniform, neutralizing background of electrons and, at densities 
above neutron drip, in a sea of neutrons. We assume that the neutron and proton 
number densities n n and n p are flat outside and inside the drop, whereas the elec- 
tron number density n e is everywhere constant. (We ignore here the thickness of 
the surface layer of the drop.) In the Wigner-Seitz approximation, the electrostatic 
energy of the Wigner-Seitz cell of radius R c reads a 

E C = 4°) (l - V/3 + l u ) t (3) 

where u = (R/R c ) 3 is the volume fraction of the drop. Equilibrium with respect to 
R with fixed number densities rij (i = n,p, e) leads to 

2E C = Ef) . (4) 

By using conditions (||) and (|j) up to C^u 1 / 3 ), we obtain the fission-instability 
criterion, u = 1/8, derived by Pethick and Ravenhall.cP Herewe have noted that 
corrections to condition (0) due to the lattice are of order u.c3 Since u generally 
increases with increasing density, we can uniquely determine the density at which 
the drop becomes unstable with respect to the deformation (|l|). 

We next derive a fission-instability condition appropriate for a nuclear rod, com- 
posed of uniformly distributed neutrons and protons (proton charge density p) as 
well as having a circular section (sectional radius R) and an infinitely long axis. 
For this rod, we consider a small sectional deformation of the quadrupole type that 
is uniform in the direction of the axis. This deformation is characterized by the 
distance, measured on a horizontal section, from the axis to a point located on the 
surface with a given angle 9: 

X(6) = R(l + a + a 2 cos26 + ---) . (5) 

Here, «o is the fractional change in the mean sectional radius of the rod, and a 2 cos 29 
represents the degree of the quadrupole- type deformation. Subject to the condition 
of constant sectional area, we obtain ao = —(l/4)a^. 

Proceeding with the argument for the rod in a vacuum, we are inevitably led 
to an infinite value of the electrostatic energy per unit length. For this reason, we 
include the effect of a triangular lattice formed by the rods immersed in a uniform 
neutralizing background of electrons and in a sea of neutrons. We here again utilize 
the Wigner-Seitz approximation. The electrostatic energy for a Wigner-Seitz cell of 
sectional radius R c can then be expressed per unit length as 

E c = ^Pp^(-lnu-l + u) , (6) 

where u = (R/R c ) 2 is the volume fraction of the rod. The deformation (|5|) produces 
a change in the electrostatic energy from this Eq given by 

5E c = (irpR 2 ) 2 (-± + u\a 2 2 . (7) 
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Next, we can obtain the surface energy per unit length of the undeformed rod and 
its change due to the deformation @ as 

Ej® = 2<irE suri R , (8) 

8E S = . (9) 

By combining the instability condition 5Ec + 6E S = 0, determined by Eqs. ([?]) 
and @ , with the size-equilibrium condition (||) , in which Eq and E^ are now given 
by Eqs. (|6|) and (||), we conclude that the rod is stable with respect to the deformation 
(|5|) for < u < 1. Recall that this deformation is uniform in the direction of the axis. 
Any small deviation from this uniformity with a fixed nucleon density and volume 
of the rod, however, results in further energy loss. This is because an increase in 
the surface energy, which is proportional to the deviation, cannot be compensated 
for by a change in the Coulomb energy, which is of second order in the deviation. 
We remark in passing that deformations of the liquid-crystal type as considered 
by Pethick and PotekhinEiP are elastic, and are not accompanied by the sectional 
distortion considered here. 

In the case of a nuclear slab of thickness 2R contained in a cell of thickness 
2R C , there is no quadrupole-type deformation corresponding to (|l|) and (||). All the 
deformations that hold fixed the nucleon density and volume of the slab have already 
been examined in the context of elastic deformations by Pethick and Potekhin;ll2P 
the slab experiences a restoring force for these deformations. 

We next turn to the condition for instability with respect to proton clustering 
in uniform nuclear matter neutralized and (3 equilibrated by electrons. This condi- 
tion was obtained by BBPt$ by expanding the energy density functional E[ni(r)] 
{i = n,p,e) of the system with respect to small density fluctuations <5nj(r) around 
the homogeneous state. The contribution of first order in 5ni vanishes due to the 
equilibrium of the unperturbed homogeneous system, while the second order contri- 
bution can be described in the spirit of the extended Thomas- Fermi model for finite 
nuclei asB 

E-E = y ^ v (q)\Sn p ( q )\ 2 , (10) 

where Eq is the ground-state energy, <5n p (q) is the Fourier transform of 5n p (r), and 
v(q) is the potential of the effective interaction between protons, as given by 

.2 



v(q)=v + (3q 2 + . (11) 



Here, 



Vo = ^p_ (dn P /dn n ) 2 
dn p dfi n /dn n 

2 

P = {Bp P + 2B np ( + B nn ( ) , (13) 

C = -f^, (14) 
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with /J, n (p) the neutron (proton) chemical potential, bnm the saturation density of 

symmetric nuclear matter, and Bij the matrix determining the gradient term in 
1/3 

E, and &tf ~ 0.3n e is the Thomas-Fermi screening length of the ultrarelativistic 
electrons, modifying the Coulomb term in E. The potential v(q) takes a minimum 
value u m in at q = Q, where 

Q 2 =(^) V2 -^F, (15) 

vmin = v Q + 2(4vre 2 /?) 1 / 2 - (3k^ F . (16) 

In the energy expansion up to second order in drii, the condition that uniform 
nuclear matter becomes unstable with respect to proton clustering reads t> m i n = 
O.cP Generally, v m - m is controlled by the bulk contribution vq, which increases as 
the proton fraction (or, equivalently, the density) increases, while the gradient and 
Coulomb terms, which tend to suppress the instability, are small positive-definite 
corrections to t>o.t$ Such a density dependence of vq ensures the presence of a 
density above (below) which the matter is stable (unstable) with respect to proton 
clustering. 

For possible proton clustering in phases with nuclear rods and nuclear slabs, we 
may assume that the density fluctuations Srii leading to proton clustering are con- 
fined within the nuclei, as it requires too much energy to produce proton clustering 
in the surrounding neutron gas. Thus, it is useful to consider the effective inter- 
action between protons in the nuclei, corresponding to Eq. (|ll|) in uniform nuclear 
matter. The bulk term vq can now be calculated by substituting the neutron and 
proton number densities inside the nuclei into Eq. (|12|). We use the condition vq = 
when estimating the density at which nuclear matter in the slab or rod becomes un- 
stable with respect to proton clustering. We remark that rigorous extension of the 
proton effective interaction in uniform nuclear matter to that in the nuclei requires 
equilibrium neutron and proton distributions, which are inhomogeneous throughout 
the regions containing protons (see, e.g., Ref. |6|)), and inclusion of the gradient and 
Coulomb correction terms. This Coulomb term depends explicitly on the spatial 
scales R and R c . 

We now proceed to calculate the densities at which instabilities arise with re- 
spect to deformation-induced fission of spherical nuclei and proton clustering in both 
uniform matter and nonspherical nuclei. For this purpose, we use the compressible 
liquid-drop model for nuclei that we constructed in Ref. ||). To obtain this model, 
we modified the BBP model by using fj,p = —Cin n [Eq. (4) in Ref. |8|)] and 

£ SU rf = C 2 tanh(3.5 MeV//A 0) )#^ f p [Eq. (6) in Ref. |)], where ^ 0) is the neutron 
chemical potential in the neutron gas, and E^f is the BBP-type surface tension 
[Eq. (7) in Ref. ||)]. Here, we set the parameters C\ and C 2 as C\ = 400 MeV fm 2 
and C2 = 1, corresponding to typical values of fi^ and E SUT { used in recent studies 
(see Figs. 1 and 2 in Ref. ||)). This model, allowing for a bcc lattice of spherical 
nuclei or bubbles, a two-dimensional triangular lattice of cylindrical nuclei or bub- 
bles, and a layered lattice of planar nuclei, exhibits the following first order phase 
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transitions in the ground-state neutron star matter: spherical nuclei — > cylindrical 
nuclei — > planar nuclei — ► uniform matter (with increasing density). The associated 
discontinuities in the volume fraction u and proton fraction x of the nuclear matter 
region can be seen from Fig. 1, displaying u and x as functions of the baryon density 
n&. The important point is that the transition density, ny, ~ 0.079 fm~ 3 , from the 
phase with spherical nuclei to that with cylindrical nuclei is close to the fission-like 
instability point u ~ 1/8 or, equivalently, ~ 0.061 fm -3 . This result is consistent 
with the analysis of Pethick and Ravenhall.0) 




Fig. 1. The volume fraction u and the proton 
fraction x of the nuclear matter region, cal- 
culated for the ground-state neutron star 
matter. The solid curve represents the vol- 
ume fraction, and the dashed curve repre- 
sents the proton fraction. 



Fig. 2. The minimum value u m in and the bulk 
contribution vo of the effective potential be- 
tween protons lying in the nuclear matter 
region, calculated for the phases with cylin- 
drical nuclei and planar nuclei, as well as 
for the uniform phase. The solid curve rep- 
resents the result for w m in, and the dashed 
curves represent the results for vq. 



Figure 2 shows the results for v m i n , Eq. (|l^), and vq, Eq. Jl2]) , calculated from 
the compressible liquid-drop model used here. In this calculation, we have set ?inm = 
0.17 fm- 3 and B nn = B np = B pp = 8.05 MeV fm 2 so as to reproduce the gradient 
term in model I of Oyamatsu;tl) this term is consistent with the experimental masses 
and radii of normal nuclei. From Fig. 2 we can estimate the critical density at which 
proton clustering occurs in uniform matter to be ~ 0.120 fm~ 3 . This value agrees 
well with the phase-equilibrium point, n& ~ 0.123 fm~ 3 , which can be obtained 
from Fig. 1, a feature consistent with the result of Pethick and Ravenhall.0) The 
critical wavelength 2i:Q~ l is obtained as 2ttQ~ 1 ~ 20 fm. This is comparable to 
an internuclear spacing of ~ 2R C ; the cell size R c , calculated for the ground state 
matter with the present nuclear model, is shown in Fig. 3. This suggests the eventual 
formation of a lattice of nuclei of some form. 

We can observe from Fig. 2 that for cylindrical and planar nuclei, the bulk term 
vq increases with decreasing density and hence does not become negative, in contrast 
to the case of uniform nuclear matter. This implies that nonspherical nuclei are stable 
with respect to proton clustering. The observed density dependence of vq results from 
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the facts that vq increases with increas- 
ing proton fraction x and that, as seen 
in Fig. 1, x increases with decreasing 
density. Moreover, calculations of the 
difference between the proton chemical 
potentials in liquid (nuclear matter) and 
gas (neutron matter) regions and of the 
compressibility of the system employing 
the present nuclear model and following 
a line of argument used by BBp|) show 
that there is no indication of the drip of 
protons out of nonspherical nuclei nor 
a thermodynamic instability of the sys- 
tem including nonspherical nuclei. Con- 
sequently, we can conclude that within 
the present nuclear model, no channel 
leading to the decay of nonspherical nu- 
clei is opened by the possible instabili- 
ties analyzed here. 

Finally, we consider a situation in which phases with nuclear rods and nuclear 
slabs are formed in a neutron star. Under the subsequent spin-down of the star, 
matter including nonspherical nuclei would be compressed in the equatorial region 
and decompressed in the polar region. Additionally, accretion of material from a 
companion star, if occurring, would act to compress the matter in the entire region 
of the star. During decompression, nonspherical nuclei might not separate into finite 
nuclei until a critical metastability with respect to quantum-tunnelling nucleation of 
the usual bcc phase with roughly spherical nuclei is realized, while during compres- 
sion, they might persist up to a critical density at which uniform matter nucleates 
quantum mechanically in a metastable system or the volume fraction u becomes 
sufficiently close to unity for the system to melt into uniform matter. The possible 
presence of such critical metastability would act to enlarge the stellar region con- 
taining the cylindrical and planar nuclei, as compared with that predicted using the 
equilibrium configuration. In order to estimate the size of this region, we should 
take into account dynamical aspects of nucleation processes and finite-temperature 
effects; the latter tend to reduce the critical metastability by promoting the quan- 
tum decay of nonspherical nuclei, and, as we demonstrated in Ref. |8|), by melting 
the low-dimensional lattice of nonspherical nuclei via elastic deformations of long 
wavelengths. 
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